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ABSTRACT 

Redshift space distortions caused by galaxy peculiar velocities provide a window onto the growth 
rate of large scale structure and a method for testing general relativity. We investigate through a 
comparison of N-body simulations to various extensions of perturbation theory beyond the linear 
regime, the robustness of cosmological parameter extraction, including the gravitational growth index 
7. We find that the Kaiser formula and some perturbation theory approaches bias the growth rate 
by la or more relative to the fiducial at scales as large as fc > 0.07 /i/Mpc. This bias propagates to 
estimates of the gravitational growth index as well as fi™ and the equation of state parameter and 
presents a significant challenge to modelling redshift space distortions. We also determine an accurate 
fitting function for a combination of line of sight damping and higher order angular dependence that 
allows robust modelling of the redshift space power spectrum to substantially higher k. 
Subject headings: cosmology: cosmological parameters — cosmology: large-scale structure of Universe 
— cosmology: theory 



1. INTRODUCTION 

Large scale structure surveys provide key pieces of ev- 
idence for the accelerated expansion of the universe and 
historically have made substantial contributions to cstab- 
Hsh A-Cold Dark Matter (ACDM) as the current stan- 
dard model of cosmology. Measurements of the baryon 
acoustic oscillation (BAO) scale and the shape of the 
matter power spectrum continue to be essential probes 
of cosmology. But although the amount of cosmological 
data grows, the presence of a cosmological constant A, or 
other forms of dark energy, has yet to be confirmed by ob- 
servations over much of cosmic time (from 1 < z < 1100) 
and some highly undesirable features of ACDM, such as 
the cosmic coincidence problem and a discrepancy o f 120 
orders of magnitude in the value of A ()Carrollll2001| ) , are 
still unsolved. 

A possible alternative explanation of cosmic acceler- 
ation suggests that the nature of gravity deviates from 
general relativity on large scales, thus affecting the mea- 
surement of the distance-redshift relation and the man- 
ner in which large scale structure forms. Recently, the 
redshift space distortions (RSD) seen in large scale struc- 
ture surveys have emerged as a po werful new m ethod of 
probing such a gravitational origin ()LindeHl2008D because 
they offer a means of measuring the growth rate of large 
scale structure directly. 

Currently, the most developed probes of cosmic accel- 
eration originate from distance measurements but say 
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very little about gravitational dynamics. However, a con- 
straint on the growth rate may assist in breaking a pos- 
sible degeneracy between models with the same expan- 
sion history but differing in terms of their gravitational 
physics. Eew probes are as sensitive to gravitational dy- 
namics as RSD; weak lensing measures the integrated 
growth to some redshift, but not the growth rate. Pair- 
wise velocity statistics of objects, e.g. supernovae, have 
been proposed to explore the growth rate, but this tech- 
nique will not produce competitive constraints until a 
sufficiently high number density of objects are measured 
with a future survey such as the Large Synoptic Survey 
Telescope (Bhatta charva. 201(3 ) . 

The growth rate can be measured via redshift space 
distortions to modest precision in recent and on- 
goin g large galaxy r edshif t surveys, s uch as SDSS- 
II JAbazaiian et all 120091), VVDS (iLe Eevre et al.l 



I2005h. WiggleZ (IDrinkwater et al.l |20l6|), and BOSS 
(jEisenstein et al.l 120 lit ) 



In the next 



Schlegel et al. 


200-1 


Laureiis et al. 


120091) 



that cover z 
, surv eys such 



0.1-0.9. 
BigBOSS 



;eneration, surv eys sucli a s 

SuMIRe (ISutollMol) and'Euchd 
will provide precise measurements 
over an even greater range of redshifts up to z = 2. 

In Scc.[2]wc review the physics of redshift space distor- 
tions from the galaxy velocity field. Section |3] discusses 
methods for accurately extracting the growth rate from 
the galaxy power spectrum. To compare analytic approx- 
imations from perturbation theory to the fully non-linear 
solutions, we carry out N-body simulations in SecS) We 
then investigate in Sections [SJ |6] and [71 the bias induced 
in the growth rate, due to the effects of non-linearity in 
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the density field, as a function of the maximum wavenum- 
ber kmax- These investigations are split into three classes 
of models; the Kaiser limit in Sec. [5l the quasi-linear 
Scoccimarro ansatz in Sec. |6] and the non-linear models 
given by perturbation theory in Sec. [71 We perform a 
closer investigation of the appropriate form of the damp- 
ing function and angular dependence in the redshift space 
power spectrum in Sec. IH finding an accurate, scale de- 
pendent correction factor, and in Sec. [9] we open the pa- 
rameter space to discuss the implication of systematics 
on detecting deviations from general relativity. The con- 
clusions of this paper are presented in Sec. 1101 

2. GALAXY VELOCITY FIELDS 

On large scales, the peculiar velocities of galaxies are 
dominated by the bulk flow motions induced by the gra- 
dients of gravitational potentials. These in turn arise 
from mass density fluctuations and so the velocity field 
traces the growth rate of large scale structure. The con- 
tinuity equation gives the relationship between the den- 
sity and velocity fields and in the linear density regime 
is written as: 

V-«(x) = -/<5(x), (1) 

where the 5, the density perturbation is defined in terms 
of the density p and its mean value p such that S = {p — 
p)/p ^ 1, D is the linear growth factor so 6{a) oc D{a), 
and / = din D/ din a is the linear growth rate. 

The peculiar velocities add an extra component to the 
cosmological redshift which perturbs the real space posi- 
tions of galaxies, Xr, along the line of sight: 



Xr + {1 + z)- 



H{z) 



(2) 



where Xs is the redshift space position and x is the direc- 
tion of the line of sight. The expansion rate, or Hubble 
parameter, H{a) = din a/dt and the redshift z = — 1. 
This apparent change in position produces an additional 
anisotropic component to the power spectrum or correla- 
tion function because the amount of shifting that occurs 
is dependent on the angle made with respect to the line 
of sight. 

The growth rate / can be extracted from measure- 
ments of the mass density power spectrum , in the com- 
bination f{z) ai{z) (jPercival fc Whitell2009D . where cr% is 
the rms mass fiuctuation amplitude proportional to the 
growth factor D(a). Thus, /erg dD/dh\a. Within gen- 
eral relativity, the growth is determined by the expansion 
history H(a) (assuming negligible dark energy clustering 
or interaction). To separate the effects of the expansion 
from any mo difications to the standard gravity picture, 
ILinded (I200I introduced the gravitational growth index 
7, which parameterises the true linear growth rate / as 



(3) 



for models that are matter dominated at high redshift 
and has been shown to be accurate to the subpercent 
level (even better for /crs) for a number of classes of 
cosmological models. 

The gravitational growth index 7 provides a use- 
ful extension to the cosmological model framework, 
allowing straightforward tests of growth vs expan- 
sion, and a compact method for distinguishing be- 
tween many classes of modified gravity models. For 



instance, a ACDM model has 7=0.55, with almost no 
dependence on the dark energy equation of state w, 
wherea s for a Dvali Gabadadze Porrati (DGP) gravity 
model pvali. Gabadadze fc Porr ati 2000), 7 = 0.68 and 
for many /(i?) gravity models 7 « 0.42 today. We will 
explore therefore not only the measurement of the growth 
rate /, but also its propagation into cosmological param- 
eter estimation of the matter density fim, the dark energy 
equation of state, and the gravitational growth index 7. 

Current measurements of / have be en made from 
galax y redshift s urveys from the 2dF ([Hawkins et al.l 
[2OOI and SDSS [Samushia etall (pOTl and older mea- 
surements] at low redshift; a nd at higher redshifts us- 
ing V VDS (IGuzzo et al.l[200l and WiggleZ (iBlake et al.l 
[2011f ). These have constrained /cts at particular red- 
shifts to ~ 10% at best, but the precision and redshift 
coverage are not yet sufficient to stringently test for mod- 
ifications to gravity. Nonetheless, they are exciting pre- 
cursors to the constraints that will be produced by BOSS 
and BigBOSS. As the measurements get more precise, it 
is important to ensure that the analysis methods and 
theoretical knowledge keep pace in their accuracy. 

Lindcr (2008) demonstrated that once / was extracted 
from next generation measurements of redshift space dis- 
tortions one could place significant constraints on 7 and 
theories of gravity. Propagating errors, one sees that at 
a single redshift the relation between the uncertainties 
on / and 7 is 
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(4) 



For an experiment to determine 7 say to 0.04, we would 
require a 2% measurement of / at 2; = 0.5 or a 1% mea- 
surement at z = 1 if we knew the expansion history to 
precisely follow a ACDM concordance cosmology. With- 
out perfect knowledge of rim(a), the constraint on / 
would need to be tighter to achieve the same precison 
on 7. Thus the ability to measure accurately / from a 
survey is a crucial topic to investigate. This article inves- 
tigates the impact of various systematics on the quality 
of cosmological constraints as we approach the precision 
required to produce strong tests of our understanding of 
gravity. 

Methods for realistically extracting the growth rate 
from the measured galaxy power s p ectru m have been 
studied by, e. g., [O kvmrur a fc Jind ([201 1[) in the lin- 
ear reg ime and [Percival fc White! ([20091 ): [Jennings et al.l 
([2011a ) in the quasi- linear regime. We extend their anal- 
yses by investigating more fully the effects of including 
non-linearities and their impact on a broader set of cos- 
mological parameters that are expected to influence mea- 
surements of the growth rate. We test how well various 
models for redshift space distortions perform, focusing 
on three redshifts, z = 0, 0.5, 1 with three cuts in scale; 
fcmar = 0.07/i/Mpc to test RSD on large scales well in 
the linear regime, kmax ~ 0.1/i/Mpc, at the onset of 
non-linearity and kmax = 0.2/i/Mpc, to investigate the 
non-linear regime, although it is really distinctions such 
as these that we are aiming to probe. 

3. REDSHIFT SPACE DISTORTION THEORY 

Redshift space distortions introduce an anisotropic 
component to the power spectrum, as the peculiar ve- 
locity of the galaxy projected along the line of sight adds 
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to the cosmological rcdshift, perturbing the galaxy po- 
sitions. This occurs as a bulk effect and to be able to 
reliably extract the growth rate from it we need to mea- 
sure a statistic such as the power spectrum of density 
fluctuations or the spatial correlation function, which 
quantify the degree to which objects cluster. In this 
work, we consider the power spectrum because of its 
close relationship with theory and the ease with which 
a linear power spectrum m ay be obtained from a Boltz- 
marni code such as CAMB (jLewis. Challinor fc Lasenbvl 
|2000[ ) . Other works have utilised th e correlation functio n 
instead to investigate RSD, su ch as lGuzzo et al.l (|2008[ ): 
iTocchini-Valentini et all (|2011[ ). 

The complete, non-linear redshift space power spec- 
trum is given by: 

^^^^^ = / (^''*"" ^'^^"'^"^ + '^^"'^l + ^^""'^^^ ' 

(5) 

as presented in IScoccimarrol (|2004() where k is the 
wavevector, r = x' — x, Auz = m(x') — u(x) and u is 
the peculiar velocity field expressed in comoving coordi- 
nates such that u = Vped H . Even though we can write 
down the full expression, it involves quantities, such as 
the non- linear density contrast ^(x), that we cannot eas- 
ily link to theoretical predictions when given a set of cos- 
mological parameters, and thus the growth rate cannot 
be extracted this way. 

In the linear density limit, the above equation re- 
duces to somethin g tractable known as the Kaiser for- 
mula (jKaisedll987[ ): 

p'{k,^i) = {b + f^?fPss{k), (6) 

where ^ is the cosine of the angle made by k with respect 
to the line of sight (r in Eq.[S|), thus making an essentially 
two-dimensional power spectrum involving radial (line 
of sight) and transverse modes. This anisotropy is in 
contrast to the undistorted linear matter power spectrum 
Pss- Note h is the linear bias factor, relating the galaxy 
density fluctuations to the mass fluctuations. Equation|6] 
is a result of several assumptions imposed on the full 
relationship between real and redshift space in Eq. [51 
which we explore in Sec. [5l 

To extend redshift space distortion theory to smaller 
scales, several quasi-linear and non-linear models have 
been proposed. The full linear model for redshift space 
distortions includes a velocity streaming term in the form 
of an extra exponential damping term, 

P^(fc,/i) =e-(^'=^^")'(fe + /Ai2)2p,,, (7) 

where = \\/(Qtt^)] [ Pl dk and Pl is the linear power 
spectrum (jFisheii ri995[ ). Within the Kaiser limit, this 
multiplication/convolution accounts for the joint density 
and velocity probability distributions along the line of 
sight in Fourier/real space assuming that the pairwise 
velocity probability distribution function (PDF) is Gaus- 
si an. 

iFisherl (|1995| ) inspired others to take a similar ap- 
proach to modell ing the more non-linear Fingers-qf-God 
(FoG) effect f see lJacksonl[l972l: I Sargent fc Turnedll"977L 
for early discussions) , and an additional small scale veloc- 
ity dispersion term ayir arises from convolving the pair- 
wise velocity dispersion profile of objects within a halo 



with the linear Kaiser or streaming model. Assuming 
a Gaussian profile for random motion of galaxies within 
virialiscd structur es results in the exponen tial damping 
term proposed bv IPeacock &: Doddsl |l994D . while using 
an expon ential profile fo r the velocity dispersion, first de- 
rived by iPeebled (I1976D and th en applied to the power 
spectrum bv iPark et al.l (|1994[) . resuhs in a Lorentzian 
damping term. Such a combination is often then taken 
to be the full quasi-linear model of redshift space distor- 
tions (IDesiacaues fc Sheth|[20Tol ) 

P^(fc,Ai) =e-(/'=^'^^)V,,,(fc,Ai)(6 + //i')'F55, (8) 
where Vyir ~ cxp[— (fctr^ir.^)^] for the Gaussian profile or 

V^ir = (l -I- k^(Tl^j.fi^) for Lorentzian damping. These 
are of course the same to first order. Note that cr^ and 
a^r are not identical in general because ay quantifies 
the bulk motion of objects, for example haloes, while 
ayir aims to model small scale motion, such as that of 
bound objects within a halo. 

However, Eqs.[7]and|8]are still linear in the sense that 
the densities and velocities are assumed to be exactly 
coherent such that the bulk flow of the dark matter traces 
out the velocity divergences exactly and Pee = f^Pss 
and Pse = —fbPss- The simplest of the models that 
accounts for the additional information contained in the 
non-linear power s pectrum of ve l ocity divergences is the 
ansatz proposed in IScoccimarrol (|2004[ ) , 

P'ik,^l) = e-(/*'-'^'^")' [b^Pssik) - 2^l'Pse + fi^Pee] , 

(9) 

where Pss, Pse and Pee are the non-linear density, 
density-velocity and velocity divergence power spectra, 
defined as Pge = {\Se*\) and Pee = ( \0\^), where 6 = 
V • V. T he velocity dispersion is given by ct,^ as defined in 
Eqn. [71 [Jennings et al.l (|2011bD demonstrated that this 
recovers the linear growth rate up to fc = 0.25 /i/Mpc to 
a precision on / of 0.64%, but it is difficult to predict 
Pss T Pse and Pee in the fully non- linear regime without 
recourse to N-body simulations. 

One of the most promising avenues for treating red- 
shift space distortions is to take a perturbative approach 
to the full transformation (Eq. [5]) and expand the term 
within angular brackets. The most widely used and sim- 
plest perturbative scheme is Standard Perturbation The- 
ory (SPT), in which the e xpansion is performed in pow- 
ers of the scale factor fsee lJain fc Bert schinger 1993, for 
more details). Such schemes are expected to model the 
non-linear effects to higher accuracy than Eq. [51 since 
that model does not constitute an exact expre ssion under 
perturbation theory as previously noted by iMatsubaral 
|2008); Taruya et al. (2010.). Some caution is required; 
no single perturbative sch eme is preferred, the y each have 
their regimes of validity ([Garlson et al.|[2009() . and they 
are limited to being useful only when the situation is 
mostly linear but requires a small correction. The de- 
marcation between linear and non-linear regimes is often 
different in redshift and real space and some perturbation 
theori es have a more limited ra nge of accuracy in redshift 
space ([Scoccimarro et aLlll999f ). Some of these problems 
are addressed by using other perturbative schemes such 
as Renormalised Perturbation Theory (RPT) and closure 
theory. 

Perturbation theory adds additional terms to the 



4 



Kwan, Lewis, Lindcr 



Table 1 

List of models considered, their free parameters and mentions in other literature. We have divided this table into three sections, 
corresponding to Sections [5] |6] and [7] as we progress down the table. Note that empirical damping means that we allow the MCMC 

process to decide the amount of damping necessary for a good fit. 



Model 



Parameters 



Equation 



Reference 



Kaiser 

Streaming 
Empirical 
Non-linear Pgg 



f,b P%k,i,) = {b + fi?f PL{k) 

f,b P^Cfc./i) = e-(-'''=f"^'')^ {b + ffj.^f PLik) 

/, b, (Ti, As above but is a free parameter 

/, b As for streaming but is replaced by non-linear Pg, 



Kaiser (19871 

Fisher (1995) 
Samushia et al. (2011 1 
Blake et arT201lT 



Scoccimarro /, b 

(with linear damping) 
Scoccimarro /, b, cr^ 

(with empirical damping) 



{k, fi) 

[b^Pss{k)-2,i''Pse+ti''Pg,^ 
As above but ay is a free parameter 




Scoccimarro 



[9] &;occimarro (2004) 

jjcmiings ct al. (2011b) 



im Heavens et al. (1998) 

^coccimarroe^^L (1999) 
Matsubara (2008 ) 

[T2l Matsubara (2008) 
[13] Taruva et al. (2010) 
[T3l Taruva et al. (2010) 



SPT 



LPT 



Taruya++ 

(with linear damping) 
Taruya++ 

(with empirical damping) 



f,b 



f, b, a-y 



{k, tM) = {b + ftl'YPL +{h+ ff^'')P!; + Pi; 



LPT 



-fcni+/(/+2)M^)CT^ X 



[PSPT + ib+ fl^^fPLkHl + fif + 2)/.2)^2] 

[b^A(k,^i,f,b)+b^Bik,flJ,b)] 
As above but ay is a free parameter 



Kaiser formula to model the non-linear contributions to 
the redshift space power spectrum. The next to linear 
order in the SPT expansion of the real space power spec- 
trum is the 1-loop expression given by: 

Pik) = Pl + Pl3 + P22, (10) 

where Pl is the linear power spectrum, and Pi3,P22 
are mode coupling terms containing extra corrections to 
third order in 6k- For redshift space, these terms are: 

P^k, m) = (fo + ffi^PL + ib + ffi^)P!3 + P22 (11) 

as der i ved by iHeayens et all ()1998D : [Scoccimarro et all 
(fT99l : iMatsubaral (|2008^ and we refer the reader to these 
works for the complete expression. Both mode coupling 
terms in Eqn. [11] now contain additional redshift space 
contributions that are sensitive to the growth rate. In 
Lagrangian Perturbation Theory (LPT), this becomes 

P^(fc,^) =e-'^'(i+/(-^+2)^')'^' [{h + f^l''fPL + {b + fii^)Pt3 

+P22 + ib + V)'PLfc'(l + /(/ + 2)^i')al] , 

(12) 

as derived by IMatsubaral (|2008[ ) . An exponential damp- 
ing term naturally arises in the derivation that accounts 
for both BAO damping and RSD smearing. Although 
some of the mode coupling terms look formidable, these 
equations are not too difficult to derive and evaluate. 
However, the expressions beyond the 1-loop terms are 
much more computationally demanding; the presence of 
mode coupling produces multidimensional integrals that 
need to be performed numerically and that can be as 
time consuming as running an entire N-body simulation. 
This is far too expensive to be done in the course of a 
parameter fitting routine during which the model power 
spectrum may need to be evaluated at each point in pa- 
rameter space. 

Another popular perturbation scheme is closure the- 
ory, in which higher order N-point statistics, such as 
the bispectrum, are defined in terms of derivatives of 



lower order N-point statistics to close the system of equa- 
tions. The closure theory ex pansion of the redsh ift space 
power spectrum derived bv iTaruva et al.l (I201CID . (here- 
after known as Taruya"'"^), produces additional terms in 
the Scoccimarro formula, and Eq. [9] is now 

P'{k,fi) = e-(AA'-")' [b^Pssik) - 2fi^Psg + fi^Pgg + 

b^A{k,fi,f,b) + b^B{k,tiJ,b)] , 

(13) 

where A{k, fi, /, b) and B{k, fi, /, b) contain the extra con- 
tributions to the power spectrum arising from the cou- 
pling of the density and velocity fields (see iTaruva et al.l 
(|2010[) for the full terms), which are rather lengthy to 
calculate even in an era where computational power is 
cheap. These three models, SPT, LPT and closure the- 
ory, comprise the main perturbative approaches to mod- 
elling the RSD power spectrum. 

4. N-BODY SIMULATIONS 

Finding a robust approach to extracting growth rate 
information requires comparison to the non-linear power 
spectra as given by simulati ons. The simu lations were 
performed using GADGET2 ([Springell 120051 ) with a box 
of length 1500 Mpc/h per side and 1024^ particles and a 
particle mass of w 2.2 x IO^-^Mq/H. We chose a reference 
ACDM cosmology with flmo = 0.25, I^a = 0.75, fib = 
0.045, = 0.97, erg = 0.8, 'h = 0.72. 

There are numerous factors that affect the ability of 
a N-body simulation to accurately reproduce the growth 
of structure but we will only outline the most significant 
effects that could impact our results. The box size was 
chosen to be large enough such that the growth of the 
modes on the scales of interest could be fed with a suf- 
ficient number of larger modes via mode coupling. The 
initial conditions for the simulati ons were obtaine d us- 
ing the Zel'dovich approximation (| ZeP dovichll 1 9 70D eval- 
uated at z = 199; this is a sufficiently high redshift to 
prevent an artificial damping of the power spectrum that 
occurs when the modes are not given enough time for 
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non-linear growth ([Heitmann et al.l 120101) . We ran 20 
simulations with the same cosmological parameters but 
each simulation is initialized by a different realization of a 
Gaussian field set by the fiducial cosmology. We also ran 
a number of simulations with different box sizes, particle 
numbers and error tolerances in the gravity calculations 
to ensure numerical convergence. Each simulation con- 
tributes three projections of a redshift distorted density 
field from the three lines of sight oriente d along each 
axis o f the box, whic h we average over as in 'Guzzo et al.l 
(|2008l ): lPercival fc W hite (2009) ; Jennings et al. (2011^ 
Although the density modes of each projection shared the 
same growth history, the distribution of particles seen 
along each axis is now different and each redshift space 
power spectrum contains new information. This gives a 
final space power spectrum that has been averaged over 
20 realisations with each simulation contributing three 
projections of the redshift space power spectrum. As an 
accuracy test of our simulations, we have compared our 
real space power spectrum averaged ove r all the realisa- 
tions t o the power spectrum emulator in [Lawrence et all 
(|2010( ) and found a dispersion of les s than 2% betwee n 
the power spectra (1% with Halofit (jSmith et al.ll2003l )) 
over the scales we consider in this article. The real space 
density and redshift space power sp ectra are rn easured 
by applying the method outlined in iJind (|2005[ ) with a 
third order mass assignment scheme (triangular shaped 
cloud, TSC) to smooth the particles on to a FFT grid 
size of 2048^ for the density power spectra. The density- 
velocity and velocity power spectra, Pqs and Peg, also 
use a TSC mass assignment scheme but with an addi- 
tional step that involves dividing by the densities to re- 
move the mass weighting imposed by binning the par- 
ticles (jScoccimarrol 12004 ). Because of the sparseness of 
the particle density in certain regions of the simulation, 
we could only achieve a maximum FFT grid resolution 
of 512^ for these power spectra. 

The cosmological parameters were de rived using the 
Metropolis algorithm of CosmoMC (|Lewis fc Bridld 
120021) to facilitate a Markov Chain Monte Carlo (MCMC) 
process to calculate the likelihood distribution of these 
parameters assuming a particular model of RSD is true. 
The errors on the redshift space power spectrum were es- 
timated using ~ 1200 Gaussian realisations (depending 
on how many were required for the covariance matrix to 
converge) using the same volume but only 256^ particles 
for speed. Applying the Zel'dovich approximation on the 
Gaussian density field then gives the velocities as: 

v{x) = aHfD(j){x), (14) 

where (f> is the potential for a particle located at x. For 
the smallest range in fc, we were able to obtain the co- 
variance matrix directly from the simulations and this 
provided a useful check against the errors given by the 
Zel'dovich approximation. We found that the resultant 
best fitting values for the growth rate were almost in- 
distinguishable for 0.03 < fc < 0.07/i/Mpc. A hn- 
ear theory power spectrum Pl is produced by CAME 
usin g the Parameterised Po st Friedmann (PPF) mod- 
ule |Fang. Hu fc Lewisll2008D to treat the perturbations 
in the fluid equations for {wojWa} models, where the 
dark energy equation of state is given by 'w{a) = wq + 
Wa{l ~ '^)- Utilizing the linear theory power spectrum. 



we have implemented SPT, LPT and Taruya^^ theory 
schemes for the real and redshift space power spectra 
within CosmoMC in order to explore the full parameter 
range that affects growth. The growth rate is modelled 
either in terms of / or 7 as an additional free parameter. 
We consider a chain to have converged when R—1 < 0.01, 
where R is defined as the ratio of the variance between 
the mean of the chains and mean of the chain variance. 

5. LIMITS OF THE KAISER FORMULA 

The Kaiser formula has already been d emonstrated to 
have a lim ited range of apph c abilit y by iJennings et al.l 
(I2011bllal): lOkuniura fc Jind (|201lD and most notably 
bv iScoccimarrol (l2004f ) but despite these shortcomings 
has been applied to a number of data sets, most of 
these at low redshift, deep in the non-linear regime. 
lOkumura fc JingI ()2011[ ): IJennings eta[\ (|2011bf ) showed 
that the Kaiser formula is unable to reproduce the growth 
rate measured from a N-body simulation on all but 
the very largest scales, especially for low mass haloes. 
We revisit the proble m using the o riginal Kaiser for- 
mula as presented in iKaiseii (jl987D but then extend 
our analysis to address its various extensions such as 
the streaming model ([Fished I1995I ). and the various 
quasi-li near FoG damping models that have bee n pro- 
posed (|PeacocF&Doddi [pp IHlkiEin The 
analysis of the Kaiser limit in this section provides a 
point of comparison for the results of the higher order 
models that we analyse in Sections |6] and [71 

Figure [1] shows the marginalised la confidence limits 
on the growth rate obtained from fits to the full redshift 
space dark matter power spectrum in {fc, /i} with 20 bins 
in fi measured from an ensemble of N-body simulations 
at z — 0, 0.5, 1. The only parameters that are allowed to 
vary are the growth rate / and the linear bias, 6, while 
all others are kept fixed at their fiducial value, with the 
exception of the empirical damping model which has 
as an additional free parameter. We have checked the 
best fitting value of the growth rate remains the same 
within Icr when the linear bias is fixed to 6 = 1 up to 
^max = 0.1/i/Mpc. However, at k^^x = 0.2/i/Mpc, we 
found that the behaviour was model dependent; some- 
times this moved the fit closer to the fiducial growth rate 
but in all circumstances, the Icr region did not lie within 
the true value anyway. 

We have quoted the mean likelihoods obtained from 
CosmoMC, rather than the best fitting points, because 
these are considered to be more robust against varia- 
tions in chain length ([Lewis fc Bridlel[200l . although the 
differences are minor. The Itr intervals are calculated 
from the ID minimum credible intervals obtained from 
the posterior by using two tailed equal likelihood limits. 
Wc consider four different models for the redshift space 
power spectrum: the Kaiser limit (Eq. [6]), the stream- 
ing model (Eq. [7]), the streaming model with non- linear 
and the empirical model which has the same form of 
exponential damping as the streaming model, but a^, is 
treated as a free parameter. These models and their free 
parameters are listed in the top third of Table [H 

Each of the four models was fitted to a redshift space 
power spectrum measured from N-body simulations with 
three cuts in fc, namely kmax = 0.07, 0.1, 0.2 /i/Mpc. For 
each of these power spectra, kmin = 0.03/i/Mpc, since 
the larger modes have been contaminated by the finite 
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Figure 1. Constraints on the growth rate using the various 
Kaiser-type models of the redshift space power spectrum evalu- 
ated at 2 = (top), z = 0.5 (middle), and 2 = 1 (bottom), with 
cutoffs at hmax = 0.07, 0. 1, 0. 2 h/Mpc. All the models are biased; 
the horizontal lines show the true value of / at each redshift. The 
fitted parameters are {/, fe}, and CTu in the case of empirical damp- 
ing. Note that some points have been offset by a small amount in 
kmaa; for clarity. 



box size. This large scale cutoff was determined by com- 
paring the real space power spectrum measured from the 
simulations to a reference power s pectrum spectrum sup - 
plicd by the CosmicEmu package (jLawrence et aljr2010( ) . 
Wc have also performed additional convergence tests in 
redshift space at an extremely high redshift of z = 19, to 
confirm that on these scales we could recover the Kaiser 
limit (Eqn[6]) with our N-body redshift space power spec- 
trum. 

The performances of these Kaiser-type models are 
unimpressive regardless of scale and redshift, even for 
kmax = 0.07/i/Mpc at z = 1 where it is commonly 
thought that structure formation might be sufficiently 
linear that the Kaiser model might be appropriate. 
While the failure of the Kai ser limit is broad l y con - 
sistent with t h e findi ngs of lOkumura &: Jing] |2011[) : 
iJennings et al.l (j2011bD . in that the growth rate is un- 
der predicted by these models, it is a little surprising 
that neither the damping term nor the inclusion of the 
non-linear power spectrum affects the best fitting value 
for the growth rate. All these variants of the Kaiser 
limit assume that Peg = f^Pss and Pse = —bfPgs, and 
the bias in these models regardless of damping suggests 
that neglecting to include the non-linearity of large scale 
motions by their simplistic relationship between density 
and velocity, rather than the smearing of power from 
FoG ef fects, is the m ost imp ortant systematic at low red- 
shifts (IScoccimarrol .2004: Je nnings et al.|[2"011a[) . But it 
is not until we examine the Scoccimarro and Taruya+^ 
models in Sections IH] and [7] that we can see the signifi- 
cance of the coupling between the density and velocity 
fields; the extra terms in the Taruya"*""*" formalism affects 
scales as large as k^ax = 0.07 /i/Mpc. 

Even though the streaming model accounts for the 
coupling between the density and velocity fields with 
the factor of cxp (— /^fc^/i^cr^j , the assumption that the 
velocity distribution is Gaussian with a scale indepen- 
dent dispersion term was shown to be of limited ap- 
plicab ility by d i rect c omparison to N-body simulations 
in Sc occimarrcl (|2004D . Neither does exchanging the lin- 
ear for the non-linear matter power spectrum in the 
Kaiser formula work; this merely shifts the power spec- 
trum in fc,// downwards to smaller fc, while leaving the 
^ dependence the same. In fact, it is the functional de- 
pendence on ^ in these RSD models that is inadequate, 
as suggested by the bias in the empirical model, which 
is allowed to adjust for as much exponential damping as 
required by the simulations. This angular dependence is 
a key point, and we explore it further in Sec. |51 

Figure shows the full 2D redshift space power spec- 
tra of the models of this section (for simplicity we do 
not show the empirical model) compared to the N-body 
power spectrum. We have shown these as a function of 
k and ^, instead of in the k±-ku plane, because we find 
it useful to fit in terms of k and ^ so these figures pro- 
vide a direct point of comparison. None of these models 
are predicting the correct behaviour at small angles (i.e. 
large /i), even on large scales (i.e. small k). Looking 
at Fig. [51 we can see that the streaming and empirical 
damping models work in an average sense via a happy 
coincidence: the linear theory matter power spectrum 
undcrpredicts the amount of damping required at large 
/i, but happens to compensate for this by predicting less 
power at ^ ~ 0, i.e. the real space power spectrum is 
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Figure 2. P''(k,fj,) for z = 0,0.5,1 calculated using the Kaiser 
limit without damping (dashed purple), the streaming model (solid 
green), and the streaming model using the non-linear matter power 
spectrum (dotted blue), compared to measured redshift space 
power spectra from N-body simulations (black). The numbers on 
the black curves indicate the contour levels in log P" {k, fi). The 
hyperbolic red contours show constant kfi = 0.05/i/Mpc. FoG ef- 
fects are manifested in the amplification of the black contours at 
large k perpendicular to the line of sight (^ Ri 0). 



smaller than it ought to be. 

The inability of the Kaiser-type formulae to predict 
an adequate amount of damping in the power spectrum 
as fi approaches unity suggests that it may be worth- 
while to consider truncating the power spectrum in kfi 
instead of just scale alone. Unfortunately, we found that 
for even quite conservative cuts such as kfi < 0.05 /i/Mpc 
and k < 0.2 /i/Mpc, no substantial improvements could 
be gained: the growth rate is still biased by more than 
la despite the increase in the error bars because of the 
weaker dependence on /. Nonetheless, we would like to 
investigate the limit at which k and kfi can be truncated 
to reproduce the streaming model with reasonable con- 
straints on /, which doubles as a consistency test. We 
found that the scales considered were prohibitively large 
(no smaller than k^ax = 0.07/i/Mpc could be allowed) 
but the best fitting values of / obtained were within la 
of the fiducial. In Fig. [21 the hyperbolic red curves show 
a contour of kfi, and we find in Sec. [5] that an accurate 
fitting function in terms of kfi can allow us to extend 
robust extraction of the growth rate / to higher k. Fur- 
thermore, in Sec.lHwe show that to truly account for FoG 
non-linear effects, we must allow for a damping term - 
or rather an amplification term - on small scales acting 
perpendicularly to the line of sight. 

6. QUASI-LINEAR MODELS 

As a next step, we investigate the Scoccimarro ansatz 
which presents a simple extension to the Kaiser formula 
beyond linear theory using the same techniques presented 
in Sec. [5l We use Eq. [5] to model the redshift space 
power spectrum and extract the growth rate. To cal- 
culate Pss, Pse and Pge, we use 1-loop SPT which was 
shown to be accurate in Paa to 1% upto fc = .08 /i/Mpc 
at z = in real space by iCarlson et al.l (|2009| ) and a few 
percent for Pgg and Paa to fc = 0.1 fe/Mp c at the same 
redshift (see Fig. 8 of ICarlson et al.|[2009l for a detailed 
comparison. Note that the accuracy is comparable with 
improved PT). 

There is a minor improvement from the Kaiser formula 
to the Scoccimarro results shown in Fig. [3l which is not 
sufficient to produce an unbiased result to la. Even on 
large scales, such as for kmax < 0.07/i/Mpc, there is 
still a small difl[erence between the Kaiser limit and the 
Scoccimarro ansatz, and both mean likelihoods are still 
skewed below the fiducial. This is rather unlikely to be 
the effect of the size of the simulation because the effect 
is model dependent, and changing the RSD model helps 
to improve the fit as we shall see in Sec. [T] 

Curiously, the Scoccimarro points in black seem to im- 
prove as the cut in k^ax loosens, particularly at z = 
where the smaller scales are weighted more strongly. The 
growth rate is being overestimated as kmax increases be- 
cause the large scales underestimate / while the small 
scales tend to overestimate it. We can see this effect at 
work in Fig. [31 in which we have also explicitly split the 
fc range into two intervals, showing the high range in red 
(the low range is given by the usual kmax = 0.1 h/Mpc 
point). If the range of our fits to larger scales has an up- 
per limit in the range 0.1 < kmax < 0.2ft,/Mpc, then we 
could coincidentally obtain a growth rate that is consis- 
tent with the fiducial value. The interpretation is slightly 
more complicated with the Scoccimarro model than in 
the Kaiser limit because there are two competing effects; 
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Figure 3. As in Fig. [T] but using the Scoccimarro ansatz with 
the real space non-Hnear power spectra evaluated using SPT. We 
consider two types of damping in these fits, exponential damping 
with cr„ predicted by linear theory and allowing cr„ to be a free, 
empirical parameter. The light red squares at kmax = 0.2h/Mpc 
are the result of fitting only over the range O.f < k < 0.2h/Mpc. 



the real space power spectrum, Pss, is underestimated 
by SPT on small scales but the real space Pse and Pgg 
power spectra in SPT diverge in comparison to the N- 
body result. The former requires a growth rate to be 
larger than the fiducial value for the model to remain 
a good fit, while the latter favours a smaller value for 
the growth rate. The inconsistencies in the Scoccimarro 
ansatz are thrown into greater relief when we consider 
what the red points are telling us: fitting over the small 
and large scales jointly is not the same as fitting the small 
and large scales independently. A more familiar setting 
for this effect is the combination of two data sets, perhaps 
from different cosmological probes, that are inconsistent 
- their combined posterior distribution need not overlap 
a similar region in parameter space as their individual 
posteriors. 

In addition, the use of empirical damping is unneces- 
sary; in terms of la confidence levels, the two models are 
indistinguishable at kmax < 0.1 h/Mpc and both models 
are equally poor on smaller scales. Neither of these two 
models are able to consistently predict the correct value 
of the growth rate and both give a value that is sub- 
stantially below the fiducial. This occurs because the 
model power spectrum does not predict enough damping 
on large scales particularly where the line of sight con- 
tribution is greatest; that is the non- linear redshift space 
power spectrum has a lower amplitude than expected as 
/.J ^ 1, but on small scales it is the reverse: FoG ef- 
fects produce too much power perpendicular to the line 
of sight. This is explored in greater detai l in Sec. [51 

Comp ared to the simulations of iJennings et al.l 
(|2011b[ ). we find two areas of departure. Although our 
simulations have the same box size as theirs, the effec- 
tive volume of our simulations obtained by averaging the 
power spectra over the full ensemble of realisations is 
substantially larger. For the complete set of simulations, 
we found weaker constraints on / and more biased fits 
for the Scoccimarro model. However, by only using half 
of the simulations, the scatter has increased such that we 
could obtain results that are consistent with the fiducial 
cosmology for the Scocci marro model. Perhaps a more 
significant difference from IJennings et al.l (|2011b[ ) is that 
we have used SPT for the real space power sp e ctra in - 
stead of the fitting function of IJennings et al] (|2011al) . 
Although this offers a convenient alternative to pertur- 
bation theory, its accuracy is limited to scales larger than 
fc « 0.1 h/Mpc because of the method used to measure 
Pgg. We elaborate on this fitting function in future work. 

If we switch to using the exact forms for Pss , Psg and 
Pgg as measured from the suite of N-body simulations, we 
can separate the effects of the two approximations that 
are involved, since we are no longer dependent on the 
accuracy of SPT for the non-linear power spectra. This 
model is represented in Fig. [4] by the green contours, 
showing the redshift space power spectrum calculated 
with the Scoccimarro formula using linear theory damp- 
ing but Pss, Pse and Pgg are measured directly from the 
N-body simulations. This is what would be produced by 
SPT if it were perfectly accurate up to fc w 0.2/i/Mpc, 
but we caution that there will be a deviation from the 
true behaviour of P^{k,^) on smaller scales because the 
resolution of the FFT grid is limited by the sparseness of 
the particles when calculating Pgg. Without the simula- 
tion results being available for every model, the results 
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Figure 4. P''(k,^) for z = 0,0.5,1 as in Fig. [2] but calculated 
using the Scoccimarro formula with linear theory damping, and 
Pee, Pes and Pgs measured from simulations (green) or 1-loop SPT 
(dashed blue) , compared to the redshift space power spectrum from 
N-body simulations (black, labelled contours). 



of the Scoccimarro ansatz could be improved over SPT 
by using a better approximation to calculate the real 
space p ower spectra such as R PT or a fitting formula 
such as I Jennings et all ()2011af ). One might be tempted 
to attribute the inadequacies of the Scoccimarro formula 
to the SPT power spectra, but this is surely not the case 
at large scales where SPT is valid to 1%. An inspection of 
Fig |4] shows that on these scales, there is a distinct differ- 
ence between the Scoccimarro ansatz with linear theory 
damping and the N-body power spectrum even at z = 1 . 
The unsatisfactory results obtained in this approach with 
SPT leads us to explore further perturbative schemes in 
the next section. 

7. LIMITS OF PERTURBATION THEORY 

We now consider more sophisticated models of red- 
shift space distortions, which are contained in the lower 
third of Table ffl These arc th e SPT model pro - 
posed bv lHeavens et alj (jlQQS l: Scoc cimarro et al.l (|1999f ) 



(Eq. [Til), the LPT model of [Matsubaral (I2008D (Eg. [12]) 
and model of iTaruva et al.l (|2010[) (Eg. [T3l) . which we 
refer to as Taruya"*""*". These models rely on first per- 
forming an expansion on the transformation from real to 
redshift space, and then the density contrast is given by 
an additional perturbative scheme in real space. To eval- 
uate the terms in Eq. [T31 we have used SPT to obtain 
Pss^Pse and Pgg for ease of comparison with the previ- 
ous Scoccimarro fits. Unlike the previous models that we 
have discussed, these models contain higher order terms 
beyond up to . The values of / obtained from us- 
ing these perturbation theory approaches are shown in 
Fig. [5l Note that we calculate both the corrections sug- 
gested by the Taruya"*""*" model and the Pss , Pse and Pee 
used in Fig [5] using SPT. 

Comparing Fig. |3|and[5] (and the shapes of the power 
spectra in Fig. and [5]), it is clear that both LPT and 
SPT models of RSD perform substantially worse than 
both the Kaiser and the Scoccimarro models. In fact, 
at z = and z = 0.5 for kmax > 0.1/i/Mpc, we were 
unable to obtain any constraints on / using these mod- 
els, as the 1(7 interval was consistent with 0. Within 
la, the values of / given by LPT and SPT are indistin- 
guishable, but the shape of the LPT power spectrum in 
Fig.[B|is somewhat more reasonable, despite the damping 
term in LPT having a stronger dependency with fi. The 
bias of the growth rate is partially due to the strength of 
the constraints; the extra information contained in the 
higher order terms means that the errors on the growth 
rate are smaller. One of the advantages of considering 
terms beyond linear theory is the additional information 
contained in those extra terms that can break the degen- 
eracy between /, b, and as- For instance, although the 
fits are more biased in SPT and LPT, P22 ^i^d Pfg con- 
tain different combinations of /, 6, and ag and some of 
these terms are proportional to {fag^al and {bagfal- 
But by far the dominant source of the discrepancy is 
that there are two levels of approximations that enter 
into the formalism: the density perturbation expansion 
and the mapping from real to redshift space. For the 
SPT redshift space power spectrum, the transformation 
from real to redshift space is expanded to third order and 
then the real spac e overdensities are p erturbed to 5^ for 
the 1-loop terms ([Heavens et al.|[l998D . The derivation 
of the LPT redshift space power spectrum proceeds in 
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Figure 5. Icr constraints on tiie growth rates derived from mod- 
elling the redshift space power spectrum with SPT, LPT and 
Taruya++ (with Hncar theory and empirical damping) models of 
the full redshift space power spectrum. We have considered the 
Taruya++ model with both linear theory and empirical damping 
terms but the SPT model does not contain any damping. Values 
of / obtained for SPT and LPT have not been shown at 2 = 0, 0.5 
for kmax = 0.2 because they were consistent with zero. 



a rather different manner, but the final expression can 
be related to SPT by expanding the linear theory damp- 
ing term that is pre dicted self consistently in the the- 
ory (jMatsubarall2008f ) . The overall effect is to reduce the 
accuracy of perturbation theory in redshift space relative 
to its efficacy in real space. There are also several coin- 
cidental occurrences that are advantageous to the Kaiser 
and Scoccimarro formulae that are absent in the LPT 
and SPT redshift space power spectra. 

The SPT and LPT models in fact miss important 
cross terms between the velocity and density fields that 
ar e included in the T aruya"'"^ expression, see Eq. 26 
of lTaruva et~aLl (|2010D . For this reason, the Taruya"*""*" 
power spectrum is able to model the redshift space power 
spectrum with a smaller bias in the growth rate, as seen 
in Fig. [5] on all scales and at all rcdshifts, but we found 
that the posteriors on b were biased beneath 6=1. It is 
worth noting, however, that all the redshift space power 
spectra considered in this section can be related to one 
another through the addition of extra terms in the per- 
turbation theory expansion. Much of the variation in 
the estimates of the redshift space power spectrum ob- 
tained from perturbation theory originates from a ju- 
dicious choice of terms to sum in the expansion of the 
density contrast. 

We emphasise that for the dashed pink curves in Fig. [B] 
(and the results presented in Fig. [5]) the SPT expansions 
of Pss, Pgg and Pee were used instead of the improved PT 
expressions an d this contributes to the deviation between 
our results and lTaruva et al.l (j2010( ) who found that their 
expression provides an unbiased description of the red- 
shift space power spectrum up to A: < 0.205 /i/Mpc at 
z ~ 1. Yet, a comparison of Fig. |4] and Fig. [5] reveals 
that it is the Scoccimarro ansatz that follows the N-body 
contours more faithfully than the Taruya^"'" model when 
using the exact non- linear expressions for Pss , Pse and 
Pee for both models. The solid pink curves in Fig.lHlthat 
correspond to the Taruya^"'' model with the real space 
Pss, Pse and Pge measured directly from N-body simu- 
lations are especially good near /i = as the Taruya++ 
expression uses an expansion in /i. None of the models 
though get the angular dependence right at lower red- 
shifts. 

8. ANGULAR DEPENDENCE AND DAMPING 

The damping term in the power spectrum arises from 
the velocities and also non-linearities. We first exam- 
ine the values of av obtained when we allow it to vary 
as an extra parameter when fitting for the growth rate. 
We allow for redshift variation between different snap- 
shots but have assumed that ay is independent of scale. 
Figure [7] shows the Icr constraints on cr^, for the Kaiser, 
Scoccimarro and Taruya'^^^ models with and without the 
linear bias (recall that the amount of damping required 
in LPT is predicted by the model) in comparison to the 
linear theory prediction shown by the black curve. The 
bottom error bars have been suppressed where they are 
consistent with = 0. 

For many of the models at most of the redshifts, cr.u is 
not helping to fit the power spectrum shape and so the 
values are consistent with = 0. Only the Scoccimarro 
model with kmax = 0.1 h/Mpc is consistent with the lin- 
ear theory predictions over all the redshifts considered to 
1(7. In fact, ay in the Taruya"*""^ and Kaiser models show 
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Figure 7. Values of (T„ determined from using empirical damping 
with the Kaiser limit (blue, centered) and the Scoccimarro formula 
(red, offset right) and the Taruya^^ model (light green, offset left) 
at ^ = 0,0.5,1 for kmax = 0.1 and 0.2/i/Mpc. For comparison, 
we have also shown the amount of damping predicted by linear 
theory (black curve). We have removed the lower error bars for 
those points whose la limit is consistent with zero. 




Figure 8. 2D isocontours of F(k,fi), using the Kaiser, Scocci- 
marro, and Taruya++ models of the redshift space power spectrum 
at 2 = using the exact N-body expressions for Pgg , Pge a-nd Pee ■ 
The lower right corner shows the prediction of linear theory, i.e. the 
standard exponential form that is used in the streaming model. 



Figure 6. P^(fe,/i) for z = 0,0.5,1 as in Fig. [2] but using 
SPT (blue), LPT (dot-dashed green) and Taruya++ (solid and 
dashed pink) treatments of the full rcdshift space power spectrum, 
and compared to N-body results (black, labelled). We plot the 
Taruya+"'" power spectrum with two different methods of evalu- 
ating the terms: dashed pink curves correspond to using SPT for 
PsS) Pse S'lid Pee and solid pink curves use the exact N-body power 
spectra. The SPT curves here do not include any damping terms 
and the Taruya++ curves use a linear theory damping term. 



a different redshift dependence compared to the fiducial, 
although it should be exactly degenerate with the lin- 
ear growth factor D{z) according to the definition below 
Eq. [71 This explains why the Taruya"'""'' power spectrum 
produces a growth rate that is closest to the fiducial cos- 
mology when ay is left as a free parameter; the functional 
dependence is genuinely distinct from the linear theory 
as we shall later see. 

If the damping term is poorly described by an expo- 
nential, then what is the correct functional form that 
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must be applied to these RSD models to reproduce th 
simulated redshift space power spectrum? Defining th 
damping term by /x) = P^bodyif^^ f^)/P^odeii'^^ 
we divide the non-linear redshift space power spectrui 
as measured in the simulations by the models that 
have considered so far. We caution that this functio 
can not be interpreted as the p airwise velocity PD F fc 
the excellent reasons outlined in lScoccimarrol (j2004[ ). Fc 
physical reasons, it is likely that F{k, fi) is composed of 
series of convolutions of different functions at each seal 
and the true pairwise velocity PDF is buried in this con 
bination. 

Nonetheless it is instructive to characterise the func 
tional form required for these analytic models to succee 
in describing the redshift space power spectrum. FigureD 
shows F{k, fi) for the Kaiser, Scoccimarro and Taruya"'' 
models. For comparison, we have also shown the lincE 
theory prediction for F{k, fi) - the standard cxponenti; 
in {fk^a^Y ~ the lower right corner of Fig. [51 A vagn 
similarity between these linear theory curves and thos 
of the Scoccimarro and Taruya^''' models only begins t 
appear for ji > 0.5 and k > 0.2 h/Mpc. Furthermore, th 
damping is actually an enhancement for much of the k- 
space, showing amplifications of 50-150%. As the RSI 
modelling improves, the maximum amplitudes of F{k, p 
are smaller, with Scoccimarro and Taruya++ models ha^ 
ing up to three times smaller deviations. This suggest 
that wc can build on their success by further modellin 
the non-linearities and coupling between the density an 
velocity fields, here through an empirical factor. 

Since we are dealing with redshift distortions, F{k, ^ 
must enter as a function of fc/x, i.e. fc^. We find that th 
curves of constant k can be described out to k^, = 0.5 b 



F{k,^,) = 



A 



2 ,,2 



1 + Bfc2^2 



(IE 



where A, B and C arc free parameters for each of th 
models. In general, the parameters are different fc 
each model, but our preliminary investigations suggcf 
that there is some sort of relationship between these va 
ues. Expansion of the denominator in kfj. would lead t 
(and hig her order) terms, reminiscen t of the quai 
tic terms in iTang. Kavo fc Takadal (|20lH) . In contras 
we find that the amplitude of F(k, ^) is non-negative o 
small scales but unity on large scales for /i = 0, as r( 
quired. The non-linear, FoG effects are responsible fc 
amplifying the RSD signal perpendicular to the line ( 
sight on small scales, and this is accurately captured b 
this form. Also noteworthy is the necessity of three fre 
parameters - the stretching and squashing effects cannc 
be simultaneously described with a single parameter as i 
linear theory damping. If we drop the extra term in C, "w 
find that F{k,fi) falls off too rapidly to describe regior 
of large kfi at small scales. The fitting form agrees wit 
the simulations at the percent level (modulo the small 
oscillations discussed below). Although the exact prove 
nance of Eq. [TS] is the topic of future work, the term i 
C arises because of extra contributions to the pairwis 
velocity PDF beyond the usual exponential distibutioi 
Futher work is also required to test this functional fori 
at other redshifts, for the moment we have restricted oui 
selves to considering z = 0. 
Figure [9] shows that in general F{k,ii) = 1 only whe 
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Figure 9. Angular dependence factor F{k,fi) required for Kaiser 
(top), Scoccimarro (middle) and Taruya++ (bottom) models to 
exactly describe the redshift space power spectrum at 2 = 0. The 
blue curves are values of F(fe,/i) as measiued directly from simu- 
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both k and ^ are small. In fact, on extremely large scales 
(fc < 0.05 /i/Mpc), F{k,fi) < 1, but it was not possi- 
ble to infer the correct form of damping in this regime 
because the range of fc/i considered is then necessarily 
small. Of greater concern, however, is that these curves 
display some oscillatory behaviour that raises the issue 
that incorrect modelling of RSD may contribute to some 
bias in the baryon acoustic feature. We have confirmed 
that when F{k, ji) is kept at a constant value of /z (Fig. |9] 
shows curves of constant k) that the oscillations vary at 
a level above the sample variance (the latter is no more 
than 2% on these scales and the oscillations are 6% above 
the background). Moreover, the shape of the oscillations 
vary with the RSD model indicating that a genuine bias 
may be present. This will be explored further in future 
work, requiring larger simulations. This will be explored 
further in future work, requiring larger simulations. 

9. TESTING GRAVITY 

We now consider extending the fitting process to in- 
clude the full set of parameters that can be reasonably 
constrained by RSD alone, such that {7, 6,a;m,wo} are 
now free parameters, but we keep the variation in w fixed 
at Wa = 0. The parameter ujm is defined in terms of 
the matter energy density such that ujm = ^mh"^ ■ We 
only use the redshift bins at z — 0.5,1 as appropriate 
for current and future surveys and limit ourselves to us- 
ing the least biased model of RSD out of each class of 
models that we have considered. These are the stream- 
ing model with non-linear matter power spectrum, the 
Scoccimarro formula with SPT power spectra and the 
Taruya"'"^ model. With the first two, we use linear the- 
ory damping; as the previous section showed, the con- 
straints in the Kaiser and Scoccimarro models are fairly 
insensitive to the value of tr^ used, but we have used an 
empirical damping term for the Taruya^"*" because this 
seems to minimise the bias in the derived values for /. 
We allow the linear bias to vary to soften the constraints 
as a best-case scenario, since we would expect some co- 
variance between the two in a real galaxy survey. We 
have also checked that this did not negatively affect the 
ability of any of these three RSD models to correctly 
model the growth rate. In all of the models considered, 
the linear bias is quite tightly constrained and is consis- 
tent with 6 = 1 at all redshifts. Furthermore, it is the 
shape of the redshift space power spectrum, rather than 
the amplitude, that most informs us about cosmology. 

Figure [10] shows the 1 and 2a posterior distributions 
in both the 7 — VLm (top row) and 7 — (bottom row) 
parameter spaces, using out to kmax = 0.1/i/Mpc. Re- 
call that 7 is the gravitational growth index and pro- 
vides a test of general relativity, which basically predicts 
7 = 0.55. 

Treating the power spectra measured 'At z — 0.5 and 
z = 1 individually, and restricting to fc < 0.1 ft./Mpc, we 
found that the joint constraints on 7 and fi^ were mostly 
unbiased; the posterior distributions of both the Scocci- 
marro ansatz and Taruya^^ model contained the fidu- 
cial values within their la regions, while for the Kaiser 
formula the fiducial values were within the la region at 
z = 1 and within the 2a limit at z = 0.5. The contours in 
the 7— Wo plane show similar characteristics, with slightly 
more bias. The marginalised posterior distributions on 
7 for both the Scoccimarro and Taruya"*"'' models are 



correct to within la, but these bounds are now generous 
enough to allow a wide swathe of cosmological models. 

The constraints on 7 from RSD substantially weaken 
when using high redshift slices, as expected from Eq.|4]as 
^^m(a) ^ 1- The width of the la contours on 7 also pre- 
vents us from drawing any interesting conclusions about 
gravity with RSD alone, although in a simplistic sense, 
the effective volume of the simulations is « 67.5 [Gpc/ft.]'^ 
if we ignore the contributions of larger scale modes. Low 
redshift surveys capable of accurately measuring 7, to- 
gether with accurate modelling of RSD, will be important 
if we wish to test gravity and cosmology. 

One of the advantages of using the gravitational growth 
index is that we expect to measure the same value of 7 
at each redshift for scale independent theories of gravity 
and we may combine observations from several redshift 
bins to strengthen our constraints. We have combined 
the power spectra in these redshift bins, neglecting any 
covariance between them. This gives the purple contours 
seen in Fig. [TOl But none of the models have benefitted 
from such treatment, despite the substantial improve- 
ment in the statistical error. This is partly because it 
would be preferable to fit each redshift slice individually 
since the sensitivity of the results to kmax is different at 
each redshift, but also indicates that some of these mod- 
els are becoming less valid as the redshift distribution 
becomes more non-linear. This is particularly true for 
the Scoccimarro ansatz, whose joint posteriors are the 
most inconsistent with the two individual distributions. 

Attempting to use power spectrum information out to 
kmax =0.2 /i/Mpc with any of these models severely ex- 
acerbates the bias, even as the statistical errors shrink. 
For example, the rightmost column of Fig. [TOlshows that 
at z = 1 where the redshift space density field might be 
expected to be reasonably linear, the posterior distribu- 
tions are far more than 2a away from their fiducial values 
regardless of the model. However, the Taruya"*"^ model 
could certainly benefit from using a different estimate of 
^'i5i5, Pse and Pgg on small scales. 

10. CONCLUSIONS 

With an combined volume of 67.5 [Gpc//i]^, we have 
tested the robustness of the Kaiser limit, Scoccimarro 
ansatz, Lagrangian and Standard perturbation theories 
and Taruya++ models of redshift space distortions, with 
a number of variations on the form of the damping term. 
Of all the linear, quasi-linear and non-linear models for 
the rcdshift space power spectrum that we have consid- 
ered, we found that the least biased of these for de- 
ter mining the growth rate / was the model proposed 
by iTaruva et al.l ()2010D based on the Taruya++ model 
if we use standard perturbation theory to evaluate the 
real space components and restrict the scales used to 
k < 0.1ft./Mpc. However, none of the models delivered 
accurate results if we extend them to fc = 0.2 h/Mpc. 

Even with sophisticated perturbation theory schemes, 
all the models that we have considered fail to describe 
correctly the angular dependence of the redshift space 
power spectrum in the quasi-linear regime and mischar- 
acterise the non-linear Finger-of-God effects. Instead, 
we present a fitting formula for the angular dependence 
factor that, in comparison to simulations, accounts accu- 
rately for both damping and FoG effects. Further work 
is needed to relate this to underlying theory and to test 
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Figure 10. 1 and 2 cr posterior distributions on 7 after extending the parameter set to include {7, fe, tuo} for the streaming model 
(leftmost column), the Scoccimarro model (second column), and Taruya^+ model (third column) for k < 0.1/i/Mpc. In each of these 
figures, the posteriors are coloured according to redshift, z = 0.5 (light green), z = I (red) and joint fits over both redshifts are in dark 
purple. Note the change in scales between the figures. The rightmost column shows the Kaiser (red), Scoccimarro (light green), and 
Taruya++ (dark purple) models using k < 0.2h/Mpc at 2: = 1. The fiducial parameters are marked with a star. 



its universality. One issue of concern is the possible in- 
teraction of the redshift space distortion modelling with 
extraction of the true baryon acoustic oscillation scale. 

The simulations provide the redshift space power spec- 
trum of the dark mat ter haloes. As work such as 
lOkumura fc Jingj ()20lH ) shows, the behavior of galaxies 
relative to dark matter can have added complications. 
We cannot fully account for the difference between what 
is measured in a galaxy survey and models of redshift 
space distortions until we also consider the problems of 
biased tracers, as well as wide angle redshift space dis- 
tortions within perturbation theory and additional astro- 
physical effects such as non-zero neutrino masses. Other 
areas of interest include using redshift space distortions 
to constrain the metric potentials, instead of the gravi- 
tational growth index, to probe the gravitational model. 

In terms of testing gravity from the growth rate, we 
found the results are not robust with respect to inaccu- 
rate redshift space distortion modeling. Bias in the con- 
fidence regions for the gravitational growth index 7 and 
other cosmological parameters can be appreciable, mak- 
ing it dangerous to use modes as large as fc = 0.1 /i/Mpc, 
even for a rather conservative set of cosmological pa- 
rameters. In fact, the set {/, 6, r2„i, wq}, is substantially 
smaller tha n those cosm ological parameters that must be 
considered ([Simpson fc Peacock 201(3) if we are to rule 
out general relativity with any confidence. 

This might be ameliorated by combination of RSD with 
other probes, though there is also the danger that this 
will merely reduce the statistical errors. If the bias in 
the RSD results is particularly strong, as is the case with 
kmax ~ 0.2ft,/Mpc, there is the additional risk that the 
RSD datasct will simply be inconsistent with and pull the 
results of other measurements or that the resultant dis- 



tribution will be multimodal. Including the parameters 
from our angular dependence factor, or its future forms, 
may provide a robust approach to mapping growth and 
gravity with redshift space distortions. 
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